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Abstract The equation for two reggeized gluons in the vacuum channel is generalized to take 
into account the running QCD coupling constant on the basis of the bootstrap condition for 
gluon reggeization. Both the gluon trajectory as a function of momentum and the interaction 
as a function of distance grow like log log in the ultraviolet. The resulting equation depends 
on the confinement region. With a simple parametrization of its influence by an effective 
gluon mass the pomeron intercept turns out much smaller than for a fixed coupling constant 
(the BFKL pomeron). 
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1 Introduction 

Recently much effort has been spent generahzing the BFKL pomeron [ 1 ] to include the 
running QCD coupling constant [ 2 ]. As in the first attempts by Lipatov et al. [ 3 ] the idea 
is to introduce the running constant directly into the equation for the vacuum channel on the 
basis of more or less plausible assumptions about the momentum or distance scale at which 
the gluons interact. Such an approach ignores the origin of the BFKL equation, which is, in 
fact, a direct consequence of the assumption that the gluon reggeizes. This assumption can 
be formalized in the so-called bootstrap condition [ 3 ]: in the gluon channel the two-gluon 
equation should possess a solution in the form of the gluon Regge pole. One can then verify 
that with the BFKL interaction kernel and the known gluon Regge-trajectory this condition 
is indeed satisfied. 

The bootstrap condition is the basic assumption in deriving the BFKL pomeron, which 
may be considered as its byproduct in relation to the vacuum channel. For that reason, in 
our opinion, all attempts to change the BFKL equation in order to introduce the running 
coupling constant should start with the bootstrap equation rather than with the vacuum 
channel one. One has to change simultaneously both the gluon interaction and its trajectory 
in a manner consistent with the bootstrap requirement. The present note is devoted to this 
problem. 

We first prove that the bootstrap equation admits more general solutions than the known 
BFKL one. These involve a generalized interaction kernel, parametrized by a function of 
a single momentum variable, related to the gluon Regge trajectory via a nonlinear integral 
equation. We determine both the kernel and the trajectory asymptotically from the known 
behaviour of the gluon density with the running coupling constant. As a result the trajectory 
turns out rising like Inlng^ at high momenta rather than like Ing^ as with a fixed coupling 
constant. The behaviour of both the kernel and the trajectory at low momenta (f- < 
remains undetermined and governed by the confinement. However the bootstrap equation 
guarantees that all the good properties of the interaction kernel for reggeized gluons remain 
the same as for the fixed coupling constant BFKL case: the kernel results insensitive to the 
infrared region. This leaves some hope that the results will depend only weakly on the details 
of the confinement mechanism. 

Parametrizing the confinement effects by a single parameter (the "gluon mass" m) and 
interpolating both the trajectory and the kernel between the ultraviolet and infrared regions 
by their asymptotic expressions at high momenta appropriately modified to match the low 
momentum behaviour we investigate the resulting equation for the pomeron at q = 0. Con- 
trary to the fixed constant case the pomeron singularity now becomes a pole in the j plane 
whose position depends on the single dimensionless parameter m/A > 1. In the limit m — A 
the equation becomes singular and effectively coincides with the BFKL one. With growing m, 
as variational calculations indicate, the pole moves to the left so that the intercept lies much 
closer to unity than for the BFKL pomeron. According to these variational estimates, the 
pomeron intercept (minus one) A with the running coupling constant for reasonable values 
of m/A weakly depends on this ratio and lies in the region of 0.15 - 0.17. 

2 The bootstrap equation and its solution 

The bootstrap condition is a result of the requirement that the equation for two reggeized 
gluons have a Regge pole solution in the gluon channel. Let a{q) = 1 -|- a;(g) be the gluon 
Regge trajectory and K{q,qi,q[) be the two-gluon interaction kernel. The (homogeneous) 
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equation for two reggeized gluons in colour representation R is then 

(j - 1 -a;(gi) -u;fe))V'(9i) = >^R j {<fq[/{'^T^f)K{q,qi,(^^)^{q[) (1) 

where <? = + g2 is the total momentum, j is the angular momentum and Xr is the strength 
of the interaction for a given colour representation. For colour group SU{3) in the gluon 
channel Ag = 3/2 and in the vacuum channel Ai = 2A8 = 3. Concentrating on the gluon 
channel and requiring that the solution of (1) be of the form of a Regge pole with a trajectory 
1 + Lo{q), one arrives at the bootstrap equation 

(a;(g) -a;(9i) -a;fe))V'(9i) = Ag / id''q[/i27rf)Kiq,qi,q[)i;{q[) (2) 

It is a complicated functional-integral equation which relates lo, K and V'- Following the 
known result for the BFKL pomeron we shall seek solutions of (2) with a constant ip{qi)- 
Then the bootstrap equation simplifies to a relation between u and K: 

uj{q) - u^iqi) - '^{12) = As y {Sq[/{2T:f)K{q, gi, q[) (3) 

It is easy to find that this relation can be satisfied provided K is chosen in the form 

K{q,qi,<^^) = ir]{qi)/M)+viq2)/viq2))/viqi - q'l) - ri{q) / ^{(^^2) (4) 
where the function ■q((i) is related to the trajectory (jji^q) by a nonlinear integral equation 

uj{q) = -Xsrjiq) J id^qi/i2nf)/r]iqi)viq2) (5) 

Equations (4) and (5) allow to find the kernel K for a given gluon trajectory u>{q) or vice 
versa. The BFKL pomeron results from a particular solution of (4) and (5) with 

V{q) = 2 (6) 

where g is the fixed coupling constant. 

The form (4) garantees that if 77(g) — as g ^ the kernel K remains nonsingular at 
9i = ^2 = ^Jid goes to zero as gi,g2 0. One can observe that the leading contributions 
of the two terms in (4) cancel in these limits. One can also prove that, with (4) and (5) 
satisfied, the equation for two gluons in the vacuum channel is nonsingular at gi — gi — > 0. 
To do that use the trick of [ 3 ] and represent the trajectory (5) as 

a;(g) = -2A8r?(g) J (d'gi/(27r)2)/r/(gi)(r?(gi) + r?(g2)) (7) 

Then, remembering that Ai = 2A8, in the vacuum channel Eq. (1) can be written as 

{j - l)V(gi) = Ai I {d\[/{27Tf){{v{qi)/v{qi - q[)mq'i)/v{q'i) - ^{qi)/M) + v{qi ' q'l)) 

+ (1^2)-,7(gMgi)M9'iM92)) (8) 

With rj{0) = the singularity at gi = q[ on the right-hand side is evidently cancelled between 
the terms coming from the kernel K and the trajectory uj. This cancellation has been long 
known for the BFKL pomeron. 
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3 The running coupling constant 

According to the considerations of the preceding section to generahze for the running couphng 
constant one has to find with it either the Regge trajectory and then construct the kernel 
K using (4) and (5) or directly the kernel K with the structure (4), that is, the function 77(g), 
and then find the trajectory from (5). Wc shall use the latter approach. To determine r](q) 
we consider the vacuum channel equation (8) at q = 0, which is equivalent to studying the 
gluon distribution for virtuality qf = 

U - mQi) = 2Ai / {d%/ {27rf )Mq,)/r,{q^ - q[)mq[)/v{q[) 

- + riiqi - q'l)) - ri{q)^{q'i) Kc^i?) (9) 

Let q\ — > 00. Then the dominant (logarithmic) contribution comes only from the first term, 
from the integration region < gi- In this limit the equation becomes 

(i - l)^(gi) = 2Ai j i<fq[/{2nf)eiql - g1)^(g'i)/r/(gi) (10) 

Comparing this with the well-known behavior of the gluon distribution at large gf we find 
asymptotically 

Vviq) = 9\q)/W = ^T^/bq^ W/^^ (n) 
where ^(g) is the QCD running coupling constant, 

6 = (11 - {2/2,)Nf)/A (12) 

and A 2± 200 Mev is the standard QCD parameter. Relation (11) is valid for g^ >> A^. So 
the comparison to the usual evolution equation allows to determine r/(g) for large values of 
the argument. 

The equation for two reggeized gluons (1) involves all values of g^, smaller than A^ in- 
cluded. As we have observed in Sec. 2, the equation is stable respective to possible singulari- 
ties of l/ri{q) as g ^ 0. However it does not mean that the equation does not depend on small 
values of g^ at all. At most we may expect that this dependence is weak. Thus we are faced 
with the problem of prolongation of the asymptotical expression for r/(g) (and the kernel K) 
into the region of g^ comparable to A^ and even smaller than A^. This cannot be done on the 
basis of any theoretical considerations. The region of small g^ is that of confinement where, 
in fact, the equation for two gluons itself looses sense. Rather we have to invoke some initial 
values for the function r]{q) at small g^ which effectively take into account the influence of 
confinement. 

The asymptotical expression (11) blows up at g^ = A^ (the infrared pole). A simple way 
to avoid it is to introduce "the gluon mass" m changing g^ to g^ + m^. Thus we assume for 
all g2 

l/r?(g) = 27r/6(g2 + m^) ln((g2 + m^)/A'') (13) 

In this formula A is the observable QCD parameter and m is a new parameter which incorpo- 
rates the confinement effects. Evidently to avoid the infrared pole wc have to choose m > A. 
The limiting case is when m = A and the pole appears at g = 0. As mentioned due to the 
properties of the kernel the vacuum channel equation remains valid for this especially simple 
case with no more parameters than A. As calculation show, values g^ << A^ are essential in 
the vacuum equation with m = A. Inspection of (13) then shows that this limit is equivalent 
to the BFKL pomeron with a coupling constant = 27r/6. One should keep in mind, how- 
ever, than there are no a priori theoretical arguments that might support the choice m = A 
and that the values g^ << A^ essential in this limit are completely unphysical (deep inside 
the confinement region). Thus the BFKL pomeron appears as a formal but unphysical limit 
of the equation with a running coupling constant with r]{qi) given by (13). 
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4 The pomeron at g = 

At q = the equation (1) can be rewritten in a simpler form changing the function ip 
according to 

ipiq) r]iq)ipiq) (14) 
Introduce the two-gluon potential V with a kernel in the momentum space 

V{q) = l/r/((?) (15) 

Then we obtain an equation 

(1 - 2rjiqmq) = 2Ai / (d^' /i2nf)Viq - q')i>{<^) - Ai77(0)F(g) / {cf^ /{2T,f)V{q'mq') 



(16) 

It has the form of the stationary Schroedinger equation in two dimensions 

Hil^ = Ei) (17) 

where the energy is E = 1 — j and the Hamiltonian is 

H = -2u{q)-2\]V + \ir]{Q)\V ><V\ (18) 

The rightmost singularity in j corresponds to the ground state energy of the Hamiltonian 
(18). 

In the Hamiltonian (18) the Regge trajectory (jj{q) plays the role of the kinetic energy. 
The interaction consists of two parts. The first one is a normal local pair interaction with a 
potential 

V{r) = J{(fq/{2n)'^)V{q)expiqr (19) 

The second part (the last term in (18)) is a separable interaction, which appears in our case 
because according to (13) 7/(0) / unless m = A. With r]{q) given by (13) both Lo{q) and 
V are completely determined by (5) and (15). However their explicit expressions in terms of 
momentum space integrals are rather complicated and can be evaluated only by numerical 
calculations. To have some idea about the behaviour of uj{q) and V{r) we therefore start 
by studying their asymptotical expressions for small and large values of their corresponding 
arguments. 

In the following we assume A = 1. 

Let us begin with the Regge trajectory. Its explicit expression is given by (5) and (13): 

u;(g) = -(A8/27r5)(g^+m^)ln(g2+m^) j (fqi/{q'i+m?)lii(ql+m^){ql+m^)ln(ql+m^) (20) 

with qi + q2 = q- Evidently Lo{q) < for all g. At g = it is finite: 

u;{0) = -(A8/26)(l + m2lnm^E^(-lnm^)) (21) 

Note that at m = 1 the second term in (21) is zero and i^(0) = —\%/2b. The asymptotical 
behaviour at >> 1 is provided by the integration over gf << g2 or << 'f' in (20). To 
study it the form (7) is more convenient. Simple calculations give the asymptotical expression 

a;(g)g2»i ~ -(Ag/^) ln(lngV Inm^) (22) 

where we have retained the Inlnm^ term to show the divergence at m = 1. As in the BFKL 
case this divergence is cancelled by the contribution from the interaction in Eq. (1). 
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The potential V{r) is explicitly given by 

V{r) = (l/27r6) J q ex.p iqr / {q^ + rri^) ln(g2 + m^) (23) 
It goes to zero as r ^ oo 

V{r = oo) = (24) 

The asymptotical behaviour at small r is determined by the behavior of 77(g) at large g^. 
Calculations give 

V{r)r«i {I /2b) ln(ln(l/r2)/ In m^) (25) 

As mentioned, in Eq.(l) the terms proportional to Inlnm^ from (22) and (25) cancel due to 

the relation Ai = 2A8- 

As one observes, the Hamiltonian (18) possesses quite normal features unless m = 1. 
The kinetic energy —2u}{q) results positive and steadily grows with g^ from its value Eq at 
g = equal to —2uj{0) and given by (21), which only leads to a constant energy shift. The 
potential energy —2XiV{r) is negative which means attraction between the gluons, partially 
compensated by the repulsive separable term. All interaction goes to zero at large values of 
r. At r ^ the potential goes to minus infinity, although the singularity at r = is quite 
weak. With such an Hamiltonian we expect normal bound states at negative values of the 
shifted energy E = E — Eq which means poles in the j-plane. 

The obtained asymptotic expressions for a;(g) and V{r) are of course valid only for large 
values of g and 1/r respectively. The exact form of both kinetic and interaction energy has to 
be calculated numerically from (20) and (23) for a given m. We know, however, that on the 
one side, the equation (1) is not very sensitive to infrared behaviour of io{q) and V{r) and, on 
the other hand, the initial function r/(g) can only be well determined at g^ >> 1, so that only 
the asymptotic expressions for u{q) at q^ » 1 and V{r) at 1/r^ >> 1 are well established. 
For these reasons, instead of using the exact values of uj{q) and V{r) folllowing from (20) and 
(23), we shall use their asymptotical expressions (22) and (25) for all g and r respectively, 
appropriately modifying them to take into account the known values of u;(0) and V{r = 00). 
Namely we shall use Eq. (17) with E changed to E — Eq and with the expressions 

a; (g) = - ( As /b) ln(ln(g2 + ) / In ) (26) 

F(r) = (l/26)ln(ln(l/r2 + m2)/lnm2) (27) 

which agree with (22) and (25) and go to zero at g = and 1/r = 0, respectively. The 
constant term Inlnm^ is in fact cancelled between the first and second term in (18). In the 
limit m = 1 the last (separable) term in (18) disappears and the Hamiltonian takes a simple 
form 

H = (2A8/6)(lnln(l + g^) - lnln(l + l/r^)) (28) 

As mentioned, in this limit small values of g and 1/r turn out to be essential. Then the 
Hamiltonian further simplifies to 

H = (2A8/6)(lng2 + Inr^) (29) 
which is nothing but the BFKL Hamiltonian up to an energy shift. 

5 Variational calculations of the pomeron intercept 

For all its nice look, Eq. (17) with kinetic and interaction terms given by (26) and (27) 
is much more difficult to solve than the BFKL one due to the lack of scale invariance. So 
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instead of trying to solve it exactly we investigate the ground state energy by the standard 
variational technique. Evidently the lowest energy eigenvalue E = 1—j depends only on the 
dimensionless ratio of the gluon mass m to the QCD parameter A = 1. 

We have taken our trial functions as linear combinations of the two-dimensional harmonic 
oscillator functions for zero angular momentum 

N 

^(^■') = E CfcLfc(ar2) exp(-arV2) (30) 
fe=o 

with Lfc the Laguerre polinomials and a and variational parameters. Up to 10 polinomials 
have been included in the present calculations. 

The resulting ground state energies in units Xi/b for different values of m are shown in 
Table 1 (the second column). Considered as a function of m the ground state energy E{m) 
has a sharp minimum at m = 1, which can be determined exactly from the known intercept 
of the BFKL pomeron. With the growth of m the energy first goes down very steeply at 
small values of m — 1 but then flattens and beginning from m — 1 ~ 0.1 diminishes very 
slowly with m. For m greater than 2 the obtained values for the ground state energy do not 
practically depend on m up to m = 6, changing from 1.3 to 1.1. Also the optimal values 
of the parameter a are shown in Table 1 (the third column). These characterize the typical 
values of for the gluons. For m = 1 a goes to zero: as mentioned small values of are 
essential in this case, which makes it completely unphysical. However with the growth of m 
the average values of q'^ rapidly grow to become considerably greater than 1 for greater 
than 1.1. Then the picture becomes self-consistent, restricted to the region of high momenta. 
In the last column of Table 1 the values of the intercept are shown for QCD with three 
flavours (Ai/6 = 4/3). For all values of m, except the ones quite close to 1, the intercept 
results considerably smaller than for the BFKL pomeron. Its typical values are of the order 
of 0.15. It is rcmarkarblc that such values are close to the ones usually taken for the pomeron 
intercept in the semiphenomenological description of soft strong phenomena in the framework 
of the supercritical pomeron model [ 4 ]. 
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Table 





E 


a=<q^ > /A^ 


A 


1.0 


-2.04 


0. 


2.72 


1.01 


-0.518 


0.50 


0.690 


1.05 


-0.309 


1.50 


0.412 


1.1 


-0.238 


2.56 


0.317 


1.5 


-0.145 


11.3 


0.193 


2.0 


-0.129 


24.0 


0.172 


3.0 


-0.119 


54.4 


0.159 


4.0 


-0.114 


90.0 


0.152 


6.0 


-0.109 


172.0 


0.145 



Table captions 

The first column gives values of the gluon mass squared in units A^ . In the second column 
ground state pomeron energies are presented in units \\/h. The third column gives the 
optimal values of the variational parameter a which correspond to average gluon momentum 
squared in units A^. The last column gives the pomeron intercept (minus one) for the QCD 
with three flavours (6 = 9/4). 



